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Abstract 

General  block  diagram  models  are  proposed  for  PWM  DC-DC  converters  in  continuous  and 
discontinuous  conduction  modes  with  fixed  switching  frequency.  Both  current  mode  control  and 
voltage  mode  control  are  addressed  in  these  models.  Based  on  these  models,  detailed  nonlinear 
and  linearized  sampled-data  dynamics  are  derived.  Asymptotic  orbital  stability  is  analyzed. 
Audio-susceptibility  and  output  impedance  are  derived.  In  this  approach,  discontinuous  con¬ 
duction  mode  and  current  mode  control  can  be  analyzed  systematically  without  special  effort.  A 
companion  paper  [1]  addresses  the  same  issues  for  the  power  stage  of  a  PWM  DC-DC  converter. 


1  Introduction 


A  DC-DC  converter  consists  of  a  power  stage  (plant)  and  a  controller  (compensated  error  amplifier), 
as  shown  in  Fig.  1.  The  objective  of  DC-DC  conversion  is  to  convert  a  source  voltage  to  a  near¬ 
constant  output  voltage  under  disturbances  at  the  source  voltage  and  load. 


Control  signal 


Source  voltage,  Vs 


Reference  voltage,  Vr 


Output  voltage,  Vo 
Load 


Figure  1:  System  diagram  of  a  DC-DC  converter 


In  this  paper,  PWM  DC-DC  converters  are  studied.  The  typical  design  approach  of  these 
converters  is  as  follows.  The  power  stage  is  modeled  and  analyzed  by  averaging  methods  [2,  3,  4]. 
Based  on  the  averaged  model,  a  controller  is  designed.  Then  simulation  programs  are  used  to  test 
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the  closed-loop  performance. 

There  are  some  shortcomings  in  the  approach  just  mentioned.  The  averaged  models  are  ap¬ 
proximate.  The  periodic  solutions  of  PWM  converters  are  averaged  to  equilibria.  The  averaged 
models  do  not  accurately  predict  subharmonic  instability  [5],  chaotic  phenomena  [6,  7,  8,  9],  and 
steady-state  DC  offset  [10].  Moreover,  it  has  been  found  that  the  directly  obtained  averaged  models 
are  inaccurate  for  converters  operated  in  discontinuous-conduction  mode  (DCM)  or  under  current 
mode  control.  This  has  necessitated  efforts  to  obtain  more  accurate  averaged  models  for  such  cases 
[11,  12,  13,  14,  15,  16,  17,  18,  4],  Because  of  the  nonlinear  nature  of  PWM  converters,  a  finite 
number  of  simulations  does  not  necessarily  reveal  all  possible  behaviors  of  PWM  converters. 

The  controlled  switch  plays  a  central  role  in  the  converter  dynamics.  It  is  therefore  important 
to  model  the  switch  operation  accurately.  Since  the  switching  action  is  essentially  discrete,  it  is 
natural  to  consider  the  use  of  sampled-data  modeling  [19,  20,  21].  In  [19],  the  ramp  function  which 
determines  switching  action  is  approximated  by  a  constant  threshold.  In  [20],  a  general  approach  is 
proposed  without  giving  detailed  modeling  of  PWM  converters.  Extensive  numerical  calculations 
are  generally  believed  to  be  needed  for  sampled-data  modeling.  This  has  led  many  researchers  to 
make  approximations  in  their  sampled-data  models,  such  as  using  a  straight  line  approximation  for 
the  system  trajectory  within  a  sampling  period.  Although  such  approximations  can  be  useful  if  the 
switching  frequency  is  sufficiently  high,  it  has  been  shown  in  [22]  that  the  resulting  models  can  be 
inaccurate  both  qualitatively  and  quantitatively. 

This  paper  extends  previous  works  on  sampled-data  modeling  for  PWM  converters  with  fixed 
switching  frequency.  General  block  diagram  models  for  PWM  DC-DC  converters  in  continuous  and 
discontinuous  conduction  modes  are  proposed.  The  ramp  function,  which  is  crucial  in  the  dynamics, 
is  carefully  modeled.  The  models  apply  to  current  mode  control  and  voltage  mode  control,  as  well 
as  other  possible  control  schemes.  Based  on  the  block  diagram  models,  detailed  nonlinear  and 
linearized  sampled-data  sampled  dynamics  are  derived.  In  this  approach,  discontinuous  conduction 
mode  and  current  mode  control  can  be  handled  systematically,  making  the  modeling  for  these  cases 
simpler  than  with  the  use  of  averaging.  Another  important  feature  of  this  work  is  that  the  results 
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are  exact  if  the  source  and  reference  voltage  signals  are  constant.  A  companion  paper  [1]  addresses 
the  same  issues  for  the  power  stage  of  a  PWM  DC-DC  converter. 

The  performance  of  DC-DC  converter  is  evaluated  by  steady-state  analysis  and  dynamic  anal¬ 
ysis.  In  steady-state  analysis,  the  existence  and  location  of  periodic  solutions  are  of  concern.  In 
dynamic  analysis,  stability  and  transient  dynamics  can  be  studied  using  the  closed-loop  eigenvalues; 
line  regulation  and  load  regulation  can  be  studied  using  audio-susceptibility  and  output  impedance. 
Both  steady-state  and  dynamic  analysis  will  be  done  in  the  paper. 

There  exist  many  possible  operating  modes  for  the  PWM  converter  because  of  its  inherent 
nonlinear  nature.  Only  two  particular  modes  are  of  interests  in  practice:  continuous  conduction 
mode  and  discontinuous  conduction  mode.  These  two  modes  have  2  and  3  stages  respectively  in 
one  cycle,  and  the  dynamics  is  linear  within  each  stage. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  2,  a  general  block  diagram  model  for 
continuous  conduction  mode  is  proposed.  By  sampling  the  state,  input  and  output  at  the  same 
rate  as  the  switching  frequency,  a  nonlinear  sampled-data  model  can  be  derived.  The  original 
periodic  solution  is  a  fixed  point  of  the  sampled-data  model.  Existence  of  fixed  points  is  studied. 
Next,  the  linearized  sampled-data  dynamics  near  a  fixed  point  is  derived.  Local  orbital  stability  is 
then  studied;  audio-susceptibility  and  output  impedance  are  derived.  In  Sec.  3,  Similar  steps  are 
applied  to  discontinuous  conduction  mode.  In  Sec.  4,  numerous  illustrative  examples  are  given. 
Conclusions  are  given  in  Sec.  5. 

2  Continuous  Conduction  Mode  (CCM) 

2.1  Block  Diagram  Model 

In  this  subsection,  a  block  diagram  model  for  the  PWM  converter  in  CCM  is  proposed.  This  model 
is  shown  to  apply  to  current  mode  control  and  voltage  mode  control. 

The  proposed  block  diagram  model  for  the  PWM  converter  in  CCM  is  shown  in  Fig.  2.  In  the 
diagram,  A\,A2  G  JHNxN,  B\ ,  13-2  G  R.Vxl.  C,  E\ .  E<>  G  R 1  x N ,  and  D  G  R  are  constant  matrices, 
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x  G  Rv,  y  G  R  are  the  state  (of  power  stage  and  controller)  and  the  feedback  signal,  respectively. 
The  source  voltage  is  vs;  the  output  voltage  is  v0.  The  notation  vr  denotes  the  reference  signal, 
which  could  be  a  voltage  or  current  reference.  The  reference  signal  vr  is  allowed  to  be  time- varying, 
although  it  is  constant  in  most  applications.  The  signal  h(t)  is  a  T-periodic  ramp.  The  clock  has  the 
same  frequency  fs  =  1/T  as  the  ramp.  This  frequency  is  called  the  switching  frequency.  The  two 
stages  in  each  clock  period  in  CCM  are  denoted  by  Si  and  S2.  The  system  is  in  Si  immediatedly 
following  a  clock  pulse,  and  switches  to  S2  at  instants  y(t)  =  h(t). 

h{t)  =  Vl  +  (Vh-Vl)(±  mod  1) 
clock 

_  y  =  Cx  +  Du 


,  Vs  x 

U  =  (  ) 

vr 


Switch  to  Si  or  S 2 


H 

\  X 
[  v0 

=  Aix  +  B\u 
=  Eix 

S2  :  | 

\  X 

=  A2x  +  B2u 

{  Vo 

=  E2x 

Switching 

Decision 


Figure  2:  A  block  diagram  model  for  PWM  converters  in  CCM 

Assume  the  switch  and  the  diode  in  the  PWM  converter  are  ideal,  so  that  there  is  no  voltage 
drop  when  they  are  on.  Most  PWM  converters  can  be  modeled  by  Fig.  2.  The  PWM  converter 
under  current  mode  control  fits  this  model  exactly,  with  h(t)  denoting  the  slope-compensating 
ramp.  For  operation  in  voltage  mode  control,  the  system  is  switched  between  stages  when  the 
ramp  signal  h(t)  intersects  with  the  feedback  signal  y.  One  of  the  switchings  generally  has  the 
same  frequency  as  the  clock.  Therefore  the  model  in  Fig.  2  is  also  good  for  the  PWM  converter 
under  voltage  mode  control.  Typical  waveforms  in  current  and  voltage  mode  control  are  shown  in 
Fig.  3  and  Fig.  4,  respectively.  Many  other  control  schemes  (e.g.,  average  current  mode  control) 
also  fit  the  model  in  Fig.  2. 
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y(t) 


h(t) 


Si  S2  Si  S2  Si  s2 


Figure  3:  Waveforms  of  a  PWM  converter  in  CCM  under  current  mode  control 


Si  S2  Si  S2  Si  s2 


Figure  4:  Waveforms  of  a  PWM  converter  in  CCM  under  voltage  mode  control 


2.2  Nonlinear  Sampled-Data  Model 

The  switching  action  is  essentially  discrete.  To  arrive  at  a  sampled-data  model,  the  operation  of 
the  PWM  converter  within  the  ( n  +  l)-st  cycle,  t  G  [nT,  (n  +  1)T),  is  considered.  Generally  in  the 
PWM  converter,  the  switching  frequency  is  sufficiently  high  that  the  variations  in  vs  and  vr  in  a 
cycle  can  be  neglected.  Thus,  take  u  =  (vs,  vr)  to  be  constant  within  the  cycle,  and  denote  its  value 
by  un  =  ( vsn,vrn ).  Let  xn  =  x(nT)  and  von  =  v0(nT).  The  sampled-data  dynamics  which  maps 
xn  to  xn+\  is  derived  next. 

Denote  by  nT  +  dn  the  switching  instant  within  the  cycle  when  y(t)  and  h(t)  intersect.  Then 


y{nT  +  dn)  =  h(nT  +  dn) 

,  and 

Si  : 

J  X 

\  v0 

=  A\x  +  B\u 
=  E\x 

for  t  G  [nT,  nT  +  dn) 

(1) 

S2  : 

j  X 

{  Vo 

=  A2x  +  B2u 
=  E2x 

for  t  G  [nT  +  dn ,  (n  +  1)T) 

(2) 

The  two  matrices  E \  and  E2  need  not  be  the  same.  For  example,  they  can  differ  if  the  equivalent 
series  resistance  (ESR)  Rc  ^  0.  When  they  differ,  the  output  voltage  is  discontinuous.  An  example 
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of  a  discontinuous  output  voltage  waveform  is  shown  in  Fig.  5.  In  most  applications,  the  output 
voltage  of  interest  is  the  maximum,  minimum,  or  average  voltage.  So  in  the  following,  E  is  used  to 
denote  either  E\ ,  E2,  or  (E±  +  E2)/2. 


a> 

O) 


ElX(nT) 


E2x(nT) 


-• - • - ► 

nT  (n+1)T  t 


Figure  5:  A  discontinuous  output  voltage  waveform 


Using  the  dynamics  (1)  and  (2)  and  the  assumption  that  u  =  (vs,vr)  is  constant  within  the 
cycle,  one  readily  obtains  the  sampled-data  dynamics  (3)- (5)  of  the  PWM  converter  in  CCM. 
The  dynamical  equation  (3)  below  is  augmented  with  the  constraint  (4),  corresponding  to  the 
intersection  condition  y(nT  +  dn)  =  h(nT  +  dn)  defining  dn. 

%n-\- 1  = 

g(xm  un)  dn^)  — 

Von 

2.3  Existence  of  T-Periodic  Solutions 

A  T-periodic  solution  for  system  in  Fig.  2  corresponds  to  a  fixed  point  for  the  sampled-data  model. 
A  periodic  solution  x°(t)  is  sampled  as  x°(0).  Let  the  fixed  point  in  the  sampled-data  dynamics 


f  (^n  ?  dn) 

rdn  r'V 

eMT-dn)(^eA1dnXn+  eAl{dn-a)dcjBlun)+  eA2{T~a)daB2un 

Jo  Jdn 


(3) 


y(nT  +  dn)  -  h(nT  +  dn) 
rd; 


f*  an 

C( eAldnxn  +  /  eAA<]n~cr'1  daBiun)  +  Dun  -  h(dn ) 

Jo 


Exr, 


(4) 

(5) 
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(3)-(5)  be  (xn,un,  dn)  =  (x°(0),  u,  d),  where  u  =  (Vs,Vry.  Then  this  fixed  point  satisfies 


z°(0)  =  f(x°(0),u,d) 

(6) 

0 

£ 

3- 

II 

O 

(7) 

where  the  functions  /(•)  and  g(-)  are  given  in  Eqs.  (3)  and  (4).  These  N  +  1  nonlinear  equations 
(Eqs.  (6)  and  (7))  in  N  +  1  unknowns  (x°(0)  and  d)  can  be  solved  by  Newton’s  method  [19].  After 
obtaining  x°(0)  and  d ,  a  periodic  solution  x°(t)  is  obtained: 


*"(«) 


eAltx°(0)  +  /q  eAl^~a^ daB\u  for  t  G  [0,  d) 

<  eA2^~d^x°(d)  +  /j  eA2^~a^daB2U  for  t  €.  [d,T) 
x°(t  mod  T)  for  t  >  T 


(8) 


A  typical  periodic  solution  x°(t)  is  shown  Fig.  6. 


Figure  6:  A  typical  periodic  solution  x°(t)  of  a  PWM  converter  in  state  space 


Next,  consider  the  existence  of  fixed  points  in  the  sampled-data  dynamics.  Assume  all  of 
the  eigenvalues  of  A\  and  A 2  are  in  the  open  left  half  of  the  complex  plane,  then  the  matrices 
I  —  eM(T-d)eAid  anc[  1  _  eA!deA2(T-d)  are  invertil)le  [23].  From  Eq.  (6),  x°(0)  can  be  expressed  as 

s°(0)  =  (/  —  e.AAT~d) eAid^-i ^eA2{T-d)  fd  eM(d-*)d(jBiU+  fT  eMT-*)d(jB2Uj 

Jo  Jd 

Similarly,  x°(d)  can  be  expressed  as  a  function  of  d,  denoted  as  Xs(d): 

Xs(d)  =  (I -eAldeA2{-T-d))-l{eAld  T  eA2{-T-'y)daB2U+  (d  eM{d-a)daBiu)  (10) 

Jd  Jo 
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The  function  Xs  ( ■ )  has  the  following  property: 


Xs(0)  =  -AflB2u  (11) 

Xa(T)  =  -A^Bm  (12) 

So  the  N  +  1  equations,  (6)  and  (7),  reduce  to  one  equation  in  one  unknown  d: 

y(d )  -  h{d)  =  CXs(d )  +  Du  -  h{d)  =  0  (13) 

The  next  result  gives  a  sufficient  condition  for  existence  of  a  fixed  point  in  the  sampled-data 
dynamics. 

Theorem  1  Assume  that  all  of  the  eigenvalues  of  A\  and  A2  are  in  the  open  left,  half  of  the  complex 
plane.  If 

(< CAflB2u  —  Du  +  h{t)))(CAflBlU  -  Du  +  h{T~ ))  <  0  (14) 

then  there  exists  a  fixed  point  (xn,un,  dn)  =  (x°(0),  u,  d)  in  the  sampled-data  dynamics  (3)-(5). 

Proof:  From  Eq.  (13),  if 

( CAf}B2u  -  Du  +  h{ti)){C Afl B\u  -Du  +  h{T~ )) 

=  (CX(0)  +  Du  -  h(0))(CXs(T)  +  Du  —  h{T~ )) 

<  0 

then  by  the  Intermediate  Value  Theorem,  there  exists  a  solution  d  satisfying  Eq.  (13).  Hence  there 
exists  a  fixed  point  (xn,un,  dn)  =  (x°(0),u,d)  in  the  sampled-data  dynamics  (3)-(5).  □ 

Note:  Whenever  there  exists  a  T-periodic  solution  for  the  system  in  Fig.  2,  it  will  be  found  as  a 
fixed  point  in  the  sampled-data  model.  However,  a  fixed  point  in  the  the  sampled-data  model  may 
not  correspond  to  a  T-periodic  solution  in  Fig.  2.  To  explain  the  last  remark,  suppose  a  fixed  point 
(xn,un,  dn)  =  (x°(0 ),u,d)  is  obtained.  Then  a  solution  x°(t)  can  be  obtained  from  Eq.  (8).  An 
example  of  signals  y°(t)  =  Cx°(t)  +  Du  and  h(t )  is  shown  in  Fig.  7.  The  solution  y°(t),  although 
periodic  with  y°(d)  =  h(d ),  cannot  be  accepted  because  two  switchings  occur  during  t  6  (0,T). 


Figure  7:  Ineligible  signals  y°(t)  and  h(t)  in  voltage  mode  control 


2.4  Linearized  Sampled-Data  Dynamics 

The  sampled-data  dynamics  (3)-(5)  is  constrained  and  nonlinear.  Assuming  that  Jj-  =  Cx°(d~ )  — 

h(d)  /  0,  the  Implicit  Function  Theorem  allows  one  to  linearize  (3)-(5)  at  the  fixed  point  ( xn ,  un,  dn )  = 
(x°(0),  u,  d)  to  obtain  an  unconstrained  linear  sampled-data  dynamical  model.  Using  a  hat  ~  to 
denote  small  perturbations  (e.g.,  xn  =  xn  —  x°(0)),  it  follows  that 

xn+i  «  $xn  +  Tun 

v  —  Ft  \^) 

uon  —  -LJ'L,n 


where 


4*  = 


df  df  dg  t  dg 

V  0  7  / 


dxr,  ddn  ddn 


dxr 


(Xn  ,Un  ,dn)  =  (x°  (0) ,U,d) 


=  e 


=  e 


A2(T-d)/j  _  ((Ai  -  A2)x°(d)  +  (.Bi  -  B2)u)C  Ald 
C(A1x°(d)  +  B1u)-h(d) 

A2(T-d)/j  _  (i°(d  )-x°(d+))C  Ald 
Cx°(d~)  -  h(d)  j 


(16) 


r  = 


df  df  dg  J  <9# 
"  V  / 


ddn  v  <9dn 


dur. 


=  eA2(T-d) 


(/ 


3AiadaBi  - 


(x„,Un,dn)  =  (x°(0),U,d) 

x°(d~)  —  x°(d+)  rd 
Cx°(d~)  -  A(d) 


(C  [d  eAiadaB1+D))+  F  ^  eA2<7daB{17) 

Jo  Jo 
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Recall  that  un  =  ( vsn,vrn )',  and  denote  T  =  [ri,r2].  Eq.  (15)  can  be  rewritten  as 


%n+ 1  ~  T  rihsrj  +  r'2'Orn 

Von  —  Fxn 


(18) 


2.5  Stability  Analysis 

In  the  PWM  converter,  the  steady-state  operating  condition  is  a  periodic  solution ,  not  an  equilib¬ 
rium  point  as  depicted  in  the  averaging  method.  The  relevant  stability  notion  is  asymptotic  orbital 
stability.  In  the  power  electronics  literature,  the  PWM  converter  is  generally  said  to  be  either 
stable  or  unstable,  without  mentioning  orbital  stability  per  se.  The  definition  of  asymptotic  orbital 
stability  is  given  as  follows. 

Definition  1  (see,  e.g.,  [24])  Denote  by  7  the  closed  orbit  generated  by  the  periodic  solution  x°(t). 
Then  x°(t )  is  asymptotically  orbitally  stable  if  there  is  a  5  such  that 

dist[x(0),  7]  <  S  =>  ^lirri  dist[x(f),7]  =  0 

where  dist[z,7]  is  defined  as  the  smallest  distance  between  the  point  z  and  any  point  on  7. 

Recall  the  following  basic  results  on  asymptotic  orbital  stability  of  the  periodic  solution  x°  ( t ) . 

Theorem  2  [24]  The  fixed  point  x° (0)  of  the  system  (3)-(5)  is  asymptotically  stable,  or  equivalently 
the  periodic  solution  x°(t)  in  the  original  continuous-time  system  of  Fig.  2  is  asymptotically  orbitally 
stable  if  all  of  the  eigenvalues  of  are  inside  the  unit  circle  of  the  complex  plane.  Moreover,  if 
x°(t)  is  asymptotically  orbitally  stable,  then  no  eigenvalues  of  $  lies  outside  the  unit  circle. 


The  next  result  follows  readily  from  the  foregoing  fact. 


Theorem  3  If  the  periodic  solution  x°(t )  is  asymptotically  orbitally  stable,  then  the  following  in¬ 
equality  holds: 


Cx°(d+)  -  h(d) 
Cx°(d~)  -  h(d) 


<;  etr[A2-Ai]d-tr[A2]T 


(19) 
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Proof:  Suppose  the  periodic  solution  x°(t )  is  asymptotically  orbitally  stable.  Then  all  the  eigenval¬ 
ues  of  $  have  magnitude  less  than  or  equal  than  1.  Since  det[<3?]  is  the  product  of  the  eigenvalues 
of  <f>,  it  follows  that 


|  det[d>]| 


det[eAldeA2{T~d)}det[I 


(x°(d  )  —  x°(d+))C 
Cx°(d~ )  -  h(d) 


_  e-tr[A2-Ai]d+tr[A2]T 


Cx°(d+ )  -  h(d) 
Cx°(d~ )  -  h{d) 


1 


The  claim  follows.  □ 

Remark:  Generally  the  switching  period  is  so  small  that  the  right  side  of  (19)  can  be  approximated 
as  1,  resulting  in  a  condition  that  resembles  a  well-known  stability  criterion  in  current  mode  control 
[4,  for  example]: 


m  2  —  mc 
m\  —  mc 


(20) 


where  m\  and  -m2  are  the  slopes  of  current  trajectories  during  the  on  and  off  stages  respectively 
using  a  linear  approximation  [2],  and  mc  is  the  slope  of  the  compensating  ramp.  Theorem  3  differs 
from  (20)  in  its  use  of  instantaneous,  rather  than  approximate,  slope. 


2.6  Comparison  with  Stability  Analysis  Using  State-Space  Averaging 

The  duty  cycle  changes  from  cycle  to  cycle.  This  variable  determines  when  the  switch  is  on  or 
off,  and  it  is  crucial  in  the  system  dynamics.  So  this  variable  needs  to  be  modeled  properly.  The 
averaging  method,  however,  treats  the  duty  cycle  as  a  continuous-time  variable. 

Next,  stability  analysis  using  the  averaging  method  is  briefly  summarized.  For  details,  the 
reader  is  referred  to  [4].  Let  the  duty  cycle  in  the  averaging  method  be  dc,  with  nominal  value  Dc. 
Also  let 


Awe  =  AqnDc  +  Aqff(1  -  Dc) 


(21) 
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BqnDc  +  Bq ff(1  -  Dc) 


(22) 


-Save  — 


Save  =  EonDc  +  Soff(1  -  Dc 

-^ave  =  ^4-ave-^ave^s 


(23) 

(24) 


The  nominal  solution  by  the  averaging  method  is  Xave,  a  constant  instead  of  a  periodic  solution 
x°(t).  Linearized  around  the  nominal  operating  point  ( x,u,dc )  =  (Xave,  (Vs,  Vr)',  Dc),  the  system 
in  Fig.  2  has  the  following  linearized  continuous-time  dynamics 

x  ~  Aavex  +  BmeU  +  ((Aon  -  A0FF)^ave  +  (Son  -  BOFF)u)dc  (25) 

y  =  Cx  +  Du  (26) 


v0  =  Eavex  +  (Son  -  S0FF)dc 


(27) 


with 


,  _  J  Vh-Vi 
ar-  ~  ' 


(for  trailing-edge  modulation,  because  dc  =  $  Yyt ) 


(for  leading-edge  modulation,  because  dc  =  ,;/ ) 


(28) 


-  il(d)T  y-L'-'-1-  XV^GtVXXXXgj  V^*_X£>V_/  XXXVWX  UXOiUXV7XX,  ^ 

In  either  trailing-edge  modulation  (Si:  switch  on,  S2 :  off)  or  leading-edge  modulation  (Si:  off, 
S2:  on),  the  closed- loop  dynamics  can  be  further  simplified  from  Eqs.  (25)- (28)  to 


x  ~  Ax  +  Bu 

v0  =  Savex  +  (Sqn  -  EOFF)dc 


(29) 


where  (stage  Si  or  S2  can  be  either  on  or  off  stage) 


A  —  Aave  +  ((Ai  —  A2)Xave  +  (B 1  —  B‘2)u)-7 


C 


B  —  Save  +  ((Ai  —  A2)Xave  +  (Si  —  B2)u)-7 


h{d)T 
D 


h(d)T 


(30) 

(31) 


The  system  is  asymptotically  stable  if  all  of  the  eigenvalues  of  A  are  in  the  open  left  half  of  the 
complex  plane,  which  is  equivalent  to  all  of  the  eigenvalues  of  exp(AT)  being  inside  the  unit  circle. 
From  Eq.  (30),  exp(AT)  can  be  written  as 


A2(T—d)+A\d+  (~(-Al^A2)Xave+(B1-B2)u)C 
exp  [AT)  =e  Md) 


(32) 
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which  is  similar  in  form  to  $  in  Eq.  (16). 

Comparing  Eqs.  (32)  and  (16),  it  is  seen  that  the  averaging  method  agrees  with  the  sampled- 
data  method  under  the  following  conditions: 

1.  x°(d)  is  approximated  by  Xave  (true  if  the  periodic  solution  x°(t)  is  approximated  by  the 
equilibrium  Xave). 

2.  Cx°(d~)  is  much  smaller  than  h(d). 

3.  T  is  very  small,  so  that  e^x°^d  )~x°(d+))CT  ~  /  _|_  (x°(d~)  —  x°(d+))CT. 

4.  The  matrices  Ai,  A2  and  (x°(d~)  —  x°(d+))C  commute. 

Note  that  condition  4  is  unlikely  to  hold,  and  condition  2  generally  does  not  hold  for  current  mode 
control. 

2.7  Audio-Susceptibility  and  Output  Impedance 

Audio-Susceptibility  and  output  impedance  are  expressed  in  terms  of  transfer  functions  (frequency 
responses) .  They  give  information  on  the  effect  of  disturbances  at  various  frequencies  on  the  output 
voltage. 

The  audio-susceptibility  is  derived  directly  from  the  linearized  sampled-data  model,  Eq.  (18). 
It  is 

Tos{z)  =  ‘^  =  E(zI-^)~1T1  (33) 

vs\z) 

To  calculate  the  output  impedance,  add  a  fictitious  current  source  ia  (as  perturbation)  in  parallel 
with  the  load.  Then  the  state  equations  in  Eqs.  (1)  and  (2)  are  replaced  by 

51  :  x  =  A\x  +  Biivs  +  Bi2vr  +  Bi3i0  (34) 

52  :  x  =  A2x  +  B2 ivs  +  B2 2vr  +  B2zi0  (35) 

Since  iQ  is  used  as  perturbation,  the  nominal  value  of  ia  is  0.  Similar  to  the  derivation  in  Sec.  2.2, 
the  new  linearized  sampled-data  dynamics  is 

•I'll-  1  ~  T  TiDsrj  +  T'2'Vrn  T  1^3 ®on  /'QR'l 

f,  —  zp/f. 

uon  — 
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where 


r3  =  -L  (  da  )~x  da 

dion  ddn  ddn  dion  ( xn,u,dn,iorl)=(x0(0),u,d,0 ) 

=  eM^d\  [d  eAia da B13  -  X°^d  ^  ~  x°[d+^>  C  [d  eAiadaB13)  [T  ^  eMadaB23  (37) 
Jo  Cx°(d~)~h(d )  Jo  Vo  V  ' 

and  ion  is  the  sampled  perturbed  output  current.  So  the  output  impedance  is 

T00(z )  =  =  E{zl  -  $)"1r3  (38) 

Given  a  transfer  function  in  z  domain,  say  T(z),  its  effective  frequency  response  [25,  p.  93]  is 
T(e^ujT),  which  is  valid  in  the  frequency  range  |(u|  <  J. 

3  Discontinuous  Conduction  Mode  (DCM) 

3.1  Block  Diagram  Model 

There  are  three  stages  in  DCM.  The  first  two  stages  and  their  operation  are  the  same  as  in  CCM. 
The  system  is  switched  to  the  third  stage  when  the  inductor  current  reaches  zero.  Within  the 
third  stage,  il  =  0.  A  block  diagram  model  for  the  PWM  converter  in  DCM  is  shown  in  Fig.  8. 
The  matrix  F  €  RlxiV  is  chosen  such  that  Fx  =  ii-  The  remaining  notation  is  the  same  as  in 
Fig.  2. 

3.2  Nonlinear  Sampled-Data  Model 

Consider  the  operation  of  the  PWM  converter  within  the  (n  +  l)-st  cycle.  As  in  CCM,  take 
u  =  (vs,vr)  to  be  constant  within  the  cycle  and  denote  its  value  by  un  =  ( vsn,vrn ). 

Denote  by  nT  +  d\n  the  switching  instant  when  y(t)  and  h(t)  intersect  within  the  cycle.  (The 
notation  d\n,  instead  of  dgn,  is  used  for  brevity.)  Denote  by  nT  +  d2n  the  switching  instant  when 
the  inductor  current  reaches  zero.  The  three  stages  Si,  S2,  S3  within  the  cycle  are  thus 

Si  :  {  X’  =  jjX  +  B1U  fo rte[nT,nT  +  dln)  (39) 

I  »o  -  n\X 
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Switching 

Decision 


h(t)  =  Vl  +  (Vh-Vl)(±  modi) 
clock 

y  =  Cx  +  Du 


il  =  Fx 


Switch  to  S\,  Si,  or  S3 


\  X 

[  Vo 

=  Aix  +  B\u 
=  Eix 

Vs  ^ 

S2  :  | 

\  X 

=  A2x  +  B2u 

vr  y  1 

[  Vo 

=  E2x 

H 

\  X 

=  A3x  +  B3u 

{  Vo 

=  E3x 

Figure  8: 

A  block  diagram  model  for  PWM  converters  in  DCM 

s2  ■■  \ 

\  X 
{  V0 

=  A2x  +  B2u 
=  E2X 

for  t  G  [nT  +  din,  nT  +  d2n) 

S3  :  | 

\  X 
[  Vo 

=  A3x  +  B3u 
=  E3x 

for  t  G  [nT  +  d2n,  ( n  +  1)T) 

The  two  switching  conditions  are 

Cx(nT  +  din)  +  Dun  =  h(nT  +  din) 
Fx(nT  +  d2n )  =  0 


(40) 


(41) 


(42) 

(43) 


As  in  CCM,  the  output  voltage  in  DCM  can  be  discontinuous  at  the  clock  time.  Again,  E  is 
used  to  denote  Ei,  E3,  or  {E\  +  E3)/2  depending  on  which  value  of  output  voltage  is  of  interest. 
From  Eqs.  (39)-(43),  the  PWM  converter  in  DCM  has  the  following  sampled-data  dynamics: 


®?i+l  —  f(xm  un,  dn) 


—  eAi(T  —  d2n)/eA2(d2n  —  d\n)^A\dlnx  _|_ 


eMdm-<T)d(jBlUn) 
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pd2n  fT 

+  /  eA2id^-a\laB2un)  +  /  eA3^T-a)daB3un 

J din  ^  d2n 


Von  —  Fxn 


g{xn,  ^n) 


Cx{nT  +  c?in)  +  .D^n  -  /i(nT  +  din) 
Fx(nT  +  d2n ) 


(44) 

(45) 


C (eAl^lri  xn  +  /0dl”  eAlC7daBiun)  +  Dun  -  h(din) 

F^eM(d2u-dln)^eA1dlnXn  +  jdm  eAi*daBlUn )  +  jd2u-dln 


=  0 


(46) 


where  dn  =  {d\nid2rf)' .  Since  the  inductor  current  always  starts  from  0  at  the  beginning  of  a  cycle, 
another  explicit  constraint  is  Fxn  =  im  =  0,  for  any  n.  So  the  dynamics  is  ( N  —  1  (-dimensional 
instead  of  iV-dimensional. 

3.3  Linearized  Sampled-Data  Dynamics 

Linearizing  Eqs.  (44)-(46)  and  using  the  notation  o  to  denote  evaluation  at  the  fixed  point  ( xn ,  un,  dn) 
(x°(0 ),u,  (d1,d2)'),  one  has 

xn+i  «  $xn  +  T  un  ,  . 

n,  —  77  ^  V4  *  / 

uon  — 


where 


$  = 


r  = 


df 


dxn 

df 


ddn 

dg 


dd„ 


df  df  dg_  i  dg 

v  £>.7  / 


dx<n  ddn  ddn  dxn 


df  df  dg  t  dg 


dnr,  dcL  dcL  dur- 


e  A3  (T— c*2 )  gA2  {d2—d\)  gAi  di 


=  e  A3(T-d2) 


eA2(d2  rfi) (£c°(rf1  )  —  i°(d+))  X°(d2  )  —  i°(^2 


Cx°(df)  -  h(di)  0 

EeA2(*-di)^o^-)  _  3.0^+))  Fx°(d^) 


(48) 

(49) 

(50) 

(51) 

(52) 
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FeM{d2-di) 


eA3(T-d2)^eA2(d2-d1)  fdl  eA1adaBi+  rh  dleA2adaB^+  fT  *  gAzOfofy  (54) 

Jo  Jo  Jo 


CJq1  eAl<TdaB\  +  D 

F{eA^d2-d.i)  jdi  eAxod(jBi  +  jd2-d!  eA2adaB 2) 


Since  the  inductor  current  %l  starts  at  0  and  ends  at  0  in  any  cycle,  it  is  not  a  dynamic  variable. 
Thus,  the  dimension  of  the  dynamics  can  be  reduced  by  one.  Note  that  this  also  implies  that  <3? 
possesses  a  zero  eigenvalue,  i.e.  that  det [T]  =  0.  This  is  shown  in  detail  as  follows: 

det[$]  =  det[eA3{T-d2)eA2id2~dl\l  ~  [  ®°(df)  -  x°(4)  e~M{d2-dl)  {x°  (d2)  -  ®°(4))  ' 


(  d9  x-i  C  \pA1d1] 

^ddj  0  FeA2{d2-di)  >e  I 


det[eAldle^(^-^)e^(*-^)]det[/-  p^-dd  ■ 
x°{di)  —  x°{df)  e~AFd2~dl\x°(d2)  —  x°(d2))  ( 


dg  x— 1 1 


=  det[eAldleA3{T~d2)eA^d2~dl)}  ■ 


C(x°(d1  )  —  x°(di))  Ce  AFd2  dl\x°(d2  )  —  x°(d2)) 


[  ~  V“1  n 


F x°(d2  ) 


Fx°(d2)  0 

-FeAFd2~dl\x°Cdi)  -  x°{dt))  Cx°(d if)  -  h{d i) 

(Cx°(df)  -h(di))Fx°((%) 


=  det[eAldleA3{T-d2)eAFd2~dl)]det[I  -  *  *  ] 


here  *  signifies  an  irrelevant  term. 


The  stability,  audio-susceptibility  and  output  impedance  analysis  are  similar  to  the  case  in  CCM 
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and  are  omitted. 


4  Illustrative  Examples 

Example  1  ( Local  vs.  global  orbital  stability ,  [10,  p.  90])  Consider  the  boost  converter  shown  in 
Fig.  9,  where  T  =  2fis,  Vs  =  4V,  L  =  5.24 /aH,  C  =  0.2 /j,F,  R  =  16P,  k\  =  —0.1,  k 2  =  0.01, 
Vr  =  0.48F,  and  h(t)  =  ((^)  mod  1). 

Although  the  control  scheme  is  neither  voltage  nor  current  mode  control,  the  circuit  can  still 
be  expressed  in  terms  of  the  block  diagram  model  in  Fig.  2,  with  state  x  =  ( iL,vc )■ 


Ai  = 

’  0  0 

0  rU 

4-2  = 

'  0 

1 

-1 

L 

-1 

Si 

=  b2  = 

1 

L 

0 

c 

RC 

C  = 

’  —ki  -I 

1 

^2 

D  = 

0 

1 

El 

=  e2  = 

0  1 

The  periodic  solution  calculated  using  Eq.  (8)  is  shown  in  Fig.  10.  The  closed-loop  poles 
calculated  from  Eq.  (16)  are  ct(<1>)  =  0.8  ±  0.45i,  which  are  inside  the  unit  circle.  So  the  periodic 
solution  is  locally  orbitally  stable.  The  magnitude  of  the  eigenvalues  is  0.9225,  indicating  that  the 
settling  time  to  the  steady  state  may  be  long.  For  example,  let  the  initial  state  be  ( iL,vc )  =  (0.9,  8). 
The  simulated  output  voltage  is  shown  in  Fig.  11.  After  the  transient,  the  state  trajectory  goes  to 
the  periodic  solution. 

The  averaging  method  also  predicts  local  stability.  From  Sec.  2.6,  the  closed-loop  poles  predicted 
by  the  averaging  method  are 

a(A)  =  (-0.2759  ±  2.9276i)  x  105  (56) 

and  e'7^1  =  0.7887  ±  0.5230?',  close  to  the  eigenvalues  found  above  using  by  the  sampled-data 
method. 

However,  the  circuit  is  not  globally  stable  and  is  described  in  [10]  as  being  unstable  based  on 
simulation.  Since  the  averaging  method  is  for  local  (small-signal)  analysis,  it  is  not  surprising  that 
it  cannot  predict  lack  of  global  stability. 
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Figure  11:  Output  voltage  trajectory  in  Example  1  for  initial  condition  (■ il,vc )  =  (0.9,8) 

Example  2  ( Buck  converter  under  voltage  mode  phase-lead  control,  [4,  p.  346])  The  system 
diagram  is  shown  in  Fig.  12.  The  system  parameters  are  as  follows:  T  =  10 [is,  Vs  =  28V ,  R  =  3f 2, 
L  =  50 [iH,  C  =  500 fiF,  Vr  =  5V,  Gc o  =  3.7,  u>z  =  10681  rad/ sec,  up  =  91106  rad/ sec,  and 
h(t)  =  4(y  mod  1).  The  voltage  divider  gain  gvci  is  chosen  to  be  0.29465  (instead  of  1/3  in  [4])  to 
result  in  an  output  voltage  at  15E. 


Figure  12:  System  diagram  for  Example  2 


Let  the  state  x  =  ( il,vc,xc ),  where  xc  is  the  state  of  the  error  amplifier.  In  terms  of  the 
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representation  in  Fig.  2, 


Ai  =  A2  = 


0 

c 


-1 

L 

-1 

RC 


0  9vd(^p  ^  z ) 


"  1 

L 

l_ 

0 

_l 

"  0 

0 

Bi 

= 

0 

0 

b2  = 

0 

0 

0 

31 

1 

3 

0 

ooz  —  U) 

p 

C 

_  GcQUJp 

UJz 

i 

1 

o 

1 

D  = 

0 

GcQUJp 

UJZ 

Ei 

=  E2  =  \ 

0  1  0 

Solving  Eqs.  (9)  and  (13)  by  Newton’s  method  results  in  x°(0)  =  (4.3,15,-0.512)  and  d  = 
5.36  x  10~6.  The  eigenvalues  of  the  closed-loop  system  can  be  obtained  from  <t[<&]  and  they  are 
0.8096  ±  0.1154  and  0.5973.  So  the  periodic  solution  x°(t)  is  asymptotically  orbitally  stable. 

Example  3  ( Boost  converter  under  current  mode  control  with  voltage  loop  closed ,  [26])  The 
system  diagram  is  shown  in  Fig.  13,  where  fs  =  25 kHz,  Vs  =  28V,  R  =  11.211,  L  =  195 pH, 
C  =  2 mF,  R±  =  47.5 fc!2,  R2  =  2.5fcl2,  Rs  =  0. 812511,  Rf  =  72. 211,  and  Cf  =  0.23 pF.  The  source 
voltage  Vs  and  thus  the  duty  cycle  are  varied.  The  system  is  analyzed  for  two  situations:  without 
slope  compensation  [h(t)  =  0),  and  with  slope  compensation  (h(t)  =  {Rs^£T)[^  mod  1]). 

Let  the  state  x  =  {iL,vc,vcf).  The  system  matrices  in  Fig.  3.1  are 


A\ 

B i 

C 

Ei 


0 

0 

0 


0 

-1 

RCJ 

cFf 

L 


0 

0 

0 


b2  = 

— Rs 

E2  =  \ 


0 

0 


CfRi 


-Rf 

Ri 

0  1 


1 

0 


0 

0 

+ 


CfR2 


a2 


D 


0 

j_ 

c 

0 


-l 

L 

-1 

RCJ 

WRJ 


0 

0 

0 


1  _i_  22L  4_  f 
1  ^  Ri  ^  R2 


Consider  the  system  without  slope  compensation.  The  duty  cycle  is  varied  from  0.4  to  0.6  and 
ct[<F]  is  plotted  and  shown  in  Fig.  14.  One  eigenvalue  trajectory  crosses  the  unit  circle  at  Dc  =  0.498. 
The  other  two  remain  very  close  to  1.  For  Dc  >  0.498,  the  system  is  unstable. 

Next,  consider  the  system  with  slope  compensation.  The  duty  cycle  is  varied  from  0.4  to  0.7 
and  er[<F]  is  calculated.  One  eigenvalue  trajectory  crosses  the  unit  circle  at  Dc  =  0.5845.  For 
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5  Concluding  Remarks 


General  block  diagram  models  have  been  proposed  for  PWM  DC-DC  converters  in  continuous  and 
discontinuous  conduction  modes  with  fixed  switching  frequency.  Both  current  mode  control  and 
voltage  mode  control  are  addressed  by  these  models.  Based  on  the  models,  detailed  nonlinear 
and  linearized  sampled-data  dynamics  have  been  derived.  Asymptotic  orbital  stability  has  been 
analyzed,  and  audio-susceptibility  and  output  impedance  have  been  derived. 

Compared  with  the  averaging  approach,  the  sampled-data  approach  has  the  following  charac¬ 
teristics  and  advantages.  An  assumption  in  the  sampled-data  approach  is  that  the  source  voltage 
and  reference  signal  can  be  viewed  as  constant  within  each  cycle.  This  assumption  is  reasonable 
because  the  switching  frequency  is  generally  very  high.  Once  the  circuit  operations  are  understood, 
the  derivation  of  the  sampled-data  dynamics  in  various  modes  and  control  schemes  are  straight 
forward  and  elegant.  For  example,  discontinuous  conduction  mode  and  current  mode  control  can 
be  analyzed  by  this  unified  approach.  Since  fewer  assumptions  are  made  than  in  the  averaging 
approach,  the  sampled-data  approach  is  more  accurate.  The  only  numerical  extensive  procedure 
in  the  sampled-data  approach  is  to  find  the  fixed-point.  The  remaining  analysis  is  eased  by  the 
analytical  form  of  the  dynamic  models.  The  sampled-data  approach  should  be  applied  to  double 
check  closed-loop  performance  whenever  a  PWM  DC-DC  converter  is  designed. 

Some  extensions  of  this  paper  and  the  companion  paper  [1]  have  been  done  in  the  first  author’s 
dissertation  [23].  These  extensions  include  modeling  and  analysis  of  load-resonant  converters,  and 
control  designs  of  both  PWM  and  load-resonant  converters. 
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